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Abstract 

A model for a flat homogeneous and isotropic Universe composed of dark energy, dark 
matter, neutrinos, radiation and baryons is analyzed. The fields of dark matter and neutrinos 
are supposed to interact with the dark energy. The dark energy is considered to obey either 
the van der Waals or the Chaplygin equations of state. The ratio between the pressure and the 
energy density of the neutrinos varies with the red-shift simulating massive and non-relativistic 
neutrinos at small red-shifts and non-massive relativistic neutrinos at high red-shifts. The 
model can reproduce the expected red-shift behaviors of the deceleration parameter and of the 
density parameters of each constituent. 

The recent astronomical measurements of type-IA supernovae [JJ [2J [5J [TJ and the analysis of 
the power spectrum of the CMBR [5J [OJ [7J [SJ [0J provided strong evidence for a present accelerated 
expansion of the Universe [31 QJJ El Q2J Q31 03] ; the nature of the responsible entity, called dark 
energy, still remains unknown. Furthermore, the measurements of the rotation curves of spiral 
galaxies |15j as well as other astronomical experiments suggest that the luminous matter represents 
only a small amount of the massive particles of the Universe, and that the more significant amount 
is related to dark matter. That offered a new setting for cosmological models with dark energy 
and dark matter and in these contexts many interesting phenomenological models appear in the 
literature analyzing the interaction of neutrinos QH QUITS] and dark matter [1511201121112311251123] 
with dark energy. With respect to dark energy some exotic equations of state were proposed in the 
literature and among others we quote the van der Waals [25, 2G, 27] [25] [20] and the Chaplygin [50J 
I5T1 [521 155] equations of state. 

In the present work a very simple cosmological model - for a homogeneous, isotropic and flat 
Universe composed by dark matter, dark energy, baryons, radiation and neutrinos - is investigated 
where the dark energy is modeled either by the van der Waals or the Chaplygin equations of state 
and interact with neutrinos and dark matter. Units have been chosen so that 8ttG/3 = c = 1, 
whereas the metric tensor has signature (+,—,—,—). 

Let a homogeneous, isotropic and spatially flat Universe be characterized by the Robertson 
Walker metric ds 2 = dt 2 — a{t) 2 8ijdx' l dx i , where a(t) denotes the cosmic scale factor. The sources 
of the gravitational field are related to a mixture of five constituents described by the fields of dark 
energy, dark matter, baryons, neutrinos and radiation. The components of the energy-momentum 
tensor of the sources is written as 

(r"„) = diag(p, -p, -p, -p), (1) 

where p and p denote the total energy density and pressure of the sources, respectively. In terms 
of the energy densities and pressures of the constituents it follows 

P = Pb+ Pdm + Pr + Pv+ Pde, P = Pb + Pdm + Pr + Pv + Pde- (2) 
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Above the indexes (b, dm, r, v, de) refer to the baryons, dark matter, radiation, neutrinos and dark 
energy, respectively. 

The conservation law of the energy-momentum tensor T^" . v = leads to the evolution equation 
for the total energy density of the sources, namely 



p + 3-(p+p) = 0, 



(3) 



where the dot refers to a differentiation with respect to time. 

The baryons and radiation are considered as non-interacting fields so that the evolution equa- 
tions for their energy densities read 



a a 
p b + 3-p b = 0, p r + A-p r = 0, 
a a 



(4) 



once the baryons represent a pressurclcss fluid, i.e., pb = 0, and the radiation pressure is given in 
terms of its energy density by p r = p r /3. 

According to a model proposed by Wetterich [T^] the evolution equation for the energy density 
of a pressureless (pdm = 0) dark matter field which interacts with a scalar field <f> is given by 



Pdm + i-Pdm = flPdm.4 1 - 



(5) 



Here the scalar field plays the role of the dark energy and (3 is a constant which couples the fields 
of dark matter and dark energy. 

For interacting neutrinos with dark energy it is supposed that the evolution equation of the 
energy density is given by (see [17l[T8] ) 



p v + 3-(p v + p v ) = a(p v - 3p„) 
a 



(6) 



The coefficient a is connected with the mass of the neutrinos and for more details one is referred 
to [TTl [I8j and to the references therein. Here a will be consider a phenomenological coefficient 
that couples the dark energy field with the neutrinos. Note that if p v = pu/3, there is no coupling 
between the fields of dark energy and neutrinos. Moreover, it is also important to note that the 
neutrinos in the past must behave as massless particles where the relationship between the pressure 
and the energy density is p v = p v /i. Due to the coupling of the neutrinos with the scalar field 
they become massive and non-relativistic. For these reasons a barotropic equation of state for the 
neutrinos is proposed where the ratio between the pressure and the energy density w v = p v j p v , 
given in terms of the red-shift z, reads 



1_ 5 *3(1/*) _ 1 ( K 3 (l/z) 
z 2 + zK 2 {\/z) z 2 \K 2 (l/z) 



(7) 



Above K 2 (l/z) and K 3 (l/z) are modified Bessel functions of second kind. For small values of z, 
w v tends to the non-relativistic limit equal to 2/3, whereas for large values of z, w v tends to the 
relativistic limit equal to 1/3. ft is noteworthy that for red-shifts z w 10 this ratio reaches the value 
w v f=a 1/3 and the coupling between the neutrinos and the dark energy is negligible. The expression 
given in ([7|) is motivated by the equation of the specific heat of a relativistic gas (see e.g. [34]). 

The evolution equation for the energy density of the dark energy field is obtained from equations 
© through ©, yielding 



Pde 



3-(pde + Pde) 

a 



-ot(j){p v - 3p„) - f3pdm<t>- 



(8) 



The energy density and pressure of the dark energy are connected with the scalar field by <f> = 
V 1 Pde + Pde- Since the purpose of this work is to develop a phenomenological theory, it is assumed 
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Figure 1: Density parameters as functions of red-shift: van der Waals fluid (solid lines) and Chap- 
lygin fluid (dashed lines). 



that the dark energy field behaves either as a van der Waals or a Chaplygin fluid with an equation 
of state given by [2H1 [29l EQl EH E21 [33] 

_ 8w vw p vw 2 A_ 

Pvw — „ °Pvwi Pen — ; \o) 

*J Pvw Pen 

where w vw and A are positive free parameters in the van der Waals and Chaplygin equations of 
state, respectively. 

For the determination of the time evolution of the energy densities one has to close the system 
of differential equations by introducing the Friedmann equation 

dV = P . (10) 

From now on the red-shift will be used as a variable instead of time thanks to the following 
relationships 

«-:-'• s = -^< 1 + s >i <"> 

Equations ^ can be easily integrated leading to the well-known dependence of the energy 
densities of the baryons and radiation with the red-shift 

p r (z) = Pr (0)(l + z)\ Pb (z) = Pb (0)(l + zf, (12) 

whereas equations (|5|) , ([6]) and ([8]) become a system of coupled differential equations for the energy 
densities pdm, Pv and pde, namely, 

(! + z )Pdm - SPdm Rn n „v 

, ==- - ~PPdm, (13) 

ViPde +Pde)/P 

{1 + z)p' v -'i(p v +p„) 

===== = -a{p u - 3p u ), (14) 

ViPde +Pde)/P 
i l + z )p'de - ?>{Pde+Pde) „ n , , n „ x 

====== = fjp dm + a(P» - 3p„). (15) 

V(Pde +Pde)/P 

In the above equations the prime refers to a differentiation with respect to the red-shift. 

In order to solve the coupled system of differential equations (fl~3|) - fp~5|) one has to specify initial 
values for the energy densities at z = 0. The following initial values for the density parameters 
Qi(z) = pi(z)/p(z) taken from the literature (see [3S] for a review) were chosen: flde{0) = 0.72, 
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Figure 2: Density parameters as functions of red-shift for a van der Waals fluid as dark matter. 



^dm(O) = 0.229916, n 6 (0) = 5 x 1(T 2 , O r (0) = 5 x 10~ 5 , fi„(0) = 3.4 x 1(T 5 . Moreover, one has 
to specify values for the coupling parameters a and ft and for the parameters w vw and A which 
appear in the van der Waals and Chaplygin equations of state One way to fix the two last 
parameters is through the use of the value of the deceleration parameter q = 1/2 + 3p/2p at z = 0. 
Indeed, by considering q(0) = —0.55 it follows w vw = 0.33851 and A = 0.50403. For the coupling 
parameters two sets of values were chosen, namely, (a) a = 5 x 10~ 5 and ft = —5 x 10~ 5 for the 
van der Waals equation of state and (b) a = 10 _1 and ft = — 10^ 2 for the Chaplygin equation of 
state. Its is also important to note that by increasing the value of the coupling parameter a (and/ 
or ft) the transfer of energy between the dark energy and neutrinos (and/or dark matter) becomes 
more efficient. 

In Fig. 1 the density parameters are plotted as functions of the red-shift for values in the range 
< z < 10. The straight lines refer to the case where the van der Waals equation of state is used 
to describe the dark energy field whereas the dashed lines correspond to the Chaplygin equation 
of state. The two density parameters that represent the dark energy field are denoted by fl vw and 
fl c h- One can infer from this figure that the dark energy density parameter tends to zero for high 
red-shifts when the van der Waals equation of state is used, whereas it tends to a constant value 
for the Chaplygin equation of state. While for high red-shifts the van der Waals equation of state 
simulates a cosmological constant with p vw = —p vw , the pressure of the Chaplygin fluid vanishes 
indicating that it becomes another component of the dark matter field (see also the behavior of 
the pressures indicated in Fig. 4). It is also important to note that the density parameters of the 
baryons and of the dark matter increase more with the red-shift for the van der Waals equation of 
state, since there is an accentuated decrease in the density parameter of the dark energy for this 
case. Note that the density parameters of the radiation and neutrinos are very small in this range 
of the red-shift and are not represented in this figure. 

The behavior of the density parameters for the cases of the van der Waals and Chaplygin 
equations of state are shown in Figs. 2 and 3, respectively, for red-shifts in the range < z < 3000. 
One can conclude from these figures, as expected, that the density parameters of the neutrinos 
and radiation increase with the red-shift whereas those of the baryons and dark matter decrease. 
Furthermore, the equality between the "matter" and "radiation" fields occurs when z as 3000 for 
the case where the dark matter field is modeled as a van der Waals fluid and z as 4200 for the case of 
a Chaplygin fluid. This can be easily understood, since in the latter case the dark energy becomes 
dark matter for high red-shifts contributing for the density parameter of the "matter" field. 

In Fig. 4 are plotted the deceleration parameter and the ratio between the pressure and the 
energy density for both cases, the large frame corresponding to the van de Waals fluid whereas the 
small frame to the Chaplygin fluid. For both cases the deceleration parameter at z = is equal to 
q(0) = —0.55, since this value was fixed in order to find the parameters w vw and A in the equations 
of state ([9]) . The transition from the decelerated to the accelerated phase of the Universe occurs at 
zt = 0.73 and zt = 0.53 for the van der Waals and Chaplygin equations of state, respectively. It 



4 





5 



0.8 




ith interactions 
ithotit interactions 



Z 



'dm 



0.6 -\ 



0.4 



'ch 



0.2 





0, 







2 



4 



6 



8 



10 



z 



Figure 5: Density parameters as functions of red-shift for a Chaplygin fluid with and without 
interactions. 



is interesting to note that while the Chaplygin equation of state simulates a cosmological constant 
with p c h = —pch for negative red-shifts which implies an accelerated phase of the Universe in the 
future, the van der Waals equation of state leads to a positive pressure and brings the Universe to 
another decelerated phase in the future. It is noteworthy to call attention that for positive values 
of the red-shift, the solution of the coupled differential equations (fT3|) through (fl"5|) predicts that 
the van der Waals fluid behaves close to a cosmological constant with p vw « —p vw . This behavior 
docs not lead to a new transition from a decelerated to an accelerated phase in the very early 
Universe, since the energy density of the radiation field increases so that the radiation pressure 
becomes larger than that of the van der Waals fluid. For high red-shifts the Universe first becomes 
dominated by the baryon and dark matter fields and for higher red-shifts by the radiation field. 
This model does not attempt to model the inflationary period, where the inflaton field dominates 
a short rapid evolution of the Universe. 

As final remarks we call attention to the fact that one expects that the coupling between dark 
energy, dark matter and neutrinos should be weak so that the parameters a and (3 are restricted to 
small values. The difference between the parameters adopted for the van der Waals and Chaplygin 
equations of state is due to stability conditions of the non-linear coupled system of differential 
equations (fl"3|) - (fT5|) . the van der Waals equation of state being more unstable for large values 
of these parameters than the Chaplygin equation of state. In Fig. 5 we have plotted the density 
parameters as functions of the red-shift for the case where a Chaplygin equation of state is used as 
dark energy. One can infer from this figure that the decay of the dark energy density parameter 
and the increase of the dark matter density parameter with the red-shift are more pronounced when 
there exists a coupling between the fields. The density parameter of the baryons remains unchanged 
since the baryons are uncoupled. 

As final comment it is important to note that even without couplings between the fields of 
dark energy, dark matter and neutrinos, this phcnomcnological model - with the equations of state 
of van der Waals and Chaplyging as dark energy - can describe satisfactorily the evolution of a 
Universe whose constituents are dark energy, dark matter, baryons, neutrinos and radiation. 
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